Polaronic features in the optical properties of the Holstein-t-J model 
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We derive the exact solution for the optical conductivity o{uj) of one hole in the Holstein-t-J 
model in the framework of dynamical mean-field theory (DMFT) . We investigate the magnetic and 
phonon features associated with polaron formation as a function of the exchange coupling J, of the 
electron-phonon interaction A and of the temperature. Our solution directly relates the features 
of the optical conductivity to the excitations in the single-particle spectral function, revealing two 
distinct mechanisms of closing and filling of the optical pseudogap that take place upon varying the 
microscopic parameters. We show that the optical absorption at the polaron crossover is charac- 
terized by a coexistence of a magnon peak at low frequency and a broad polaronic band at higher 
frequency. An analytical expression for cr(a)) valid in the polaronic regime is presented. 



I. INTRODUCTION 

The problem of a single hole in the t-J model inter- 
acting also with the lattice degrees of freedom has re- 
cently attracted a notable interest in connection with the 
physical properties of the high-Tc cupratesJ^'^'^'^ In par- 
ticular, in parent and strongly underdoped compounds, 
angle resolved photoemission spectroscopy (ARPES) re- 
veals a low energy peak whose dispersion is well described 
by the t-J model, while its anomalously large broad- 
ness has been ascribed to incoherent multi-phonon shake- 
off processesi^i^ii A similar interplay between electron- 
electron and electron-phonon interactions should in prin- 
ciple be reflected in the optical conductivity spectra. As 
a matter of fact, the most remarkable features observed 
in the underdoped region are an ubiquitous mid-infrared 
(MIR) peak at « 0.5 eV, and a weaker peak around « 0.1 
eV, the latter being more strongly doping and tempera- 
ture dependentj^i^ Several interpretations have been pro- 
posed for the origin of these features, including midgap 
or impurity states fi^iii charge/spin stripes, ^^'^'^ polaronic 
excitationSfi^ii^'i^iiiii^ and the interaction with the an- 
tiferromagnetic background rli^ 

The first proposal of the t-J model as a suitable 
basis to discuss the optical spectra of the cuprates 
was advanced by Zhang and Ricei^i who observed 
that the 1/uj behavior— of the optical absorption 
above « 0.5 eV could be naturally associated to 
the incoherent motion in an antiferromagnetic back- 
ground. Successively, the optical conductivity of the 
t-J model has been investigated in detail using several 
techniques, such as exact-diagonalization,^^ analytical 



approximations 
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and dynamical mean-field the- 



ory (DMFT).S^i21i2ai2a However, and in spite of the above 



discussed relevance of the electron-phonon coupling, the 
optical conductivity of the t-J model in the presence of 
the lattice degrees of freedom has not been thoroughly 
investigated. Numerical calculations based on exact- 
diagonalization of finite clusters were employed for in- 
stance in Ref. 30; to evaluate cr{uj). Alternatively, the 
optical conductivity was calculated analytically in Ref. 
l3ll based on a non-crossing Born approximation, which 
is however unable to describe the polaron formation. 

In this paper we present results for the optical con- 
ductivity of a single hole in the Holstein-i- J model ob- 
tained in the framework of dynamical mean-field the- 
ory. One-hole spectral properties at zero temperature 
were discussed in a previous publication where antifer- 
romagnetic correlations were shown to enhance the ef- 
fects of the electron-phonon couplingi^ A similar result 
was found also in the antiferromagnetic phase of the 
Holstein-Hubbard model using DMFT techniquesi^ A 
serious drawback of DMFT, which is obtained as the ex- 
act solution of the lattice problem in the limit of infinite 
dimensions, is that the magnetic background is treated 
in a classical way. This, together with the fact that Trug- 
man loops are negligible in infinite dimensions, prevents 
the possibility to account for coherent hole-propagation, 
which is related to the spin-flip fluctuations. As a con- 
sequence, no Drude peak can be observed in the optical 
conductivity. On the other hand, the incoherent contri- 
butions of (j{uj) are mainly dominated by local properties, 
such as the local electron-phonon scattering and spin- 
string excitations within the magnetic polaron, which are 
well captured by this approach)^-^ 

Bearing the above limitations in mind, the aim of the 
present work is thus to focus on the incoherent part of the 
optical conductivity of a single-hole and to investigate in 
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detail its features in the different physical regimes of the 
Holstein-t- J model. The dependence of the optical spec- 
tra on the microscopic parameters is analyzed with spe- 
cial attention to the intermediate coupling region, where 
the interplay between magnetic and lattice degrees of 
freedom is strongest. We show that a crucial role is 
played by the formation of the lattice polaron, which 
drives a transfer of spectral weight towards higher fre- 
quencies, opening a pseudogap at low frequencies. Con- 
versely, starting from the polaronic phase, two different 
mechanisms can be clearly identified as being responsi- 
ble for the disappearance of the pseudogap: i) reducing 
the effective exchange energy scale suppresses the posi- 
tive feedback of magnetism on polaron formation and can 
lead to a closing of the pseudogap as the system crosses 
back to the non-polaronic regime; ii) increasing the tem- 
perature, which does not alter the lattice/magnetic in- 
terplay, leads to a filling of the pseudogap more similar 
to what is expected in purely polaronic models. In the 
immediate vicinity of the polaron crossover, the spec- 
tra are characterized by a coexistence of a magnon peak 
at low frequency and a broad polaronic band at higher 
frequency, which closely resembles the experimental sit- 
uation observed in the cuprates. 

On theoretical grounds, the definition of the optical 
conductivity of a single hole is a delicate matter which 
needs particular care. We provide an analytical deriva- 
tion which generalizes the results of Refs. [27ll33 to the 
Holstein-t-J model. This approach permits us to iden- 
tify the role of the different one-particle properties on the 
optical conductivity. Comparison with numerical data is 
also discussed, showing a good agreement between our 
findings and exact diagonalization results. 

The paper is organized as follows. In Sec. |TT]we discuss 
the exact solution in infinite dimensions for the one-hole 
Green's function of the Holstein-t- J model at finite tem- 
perature. Results for the one-particle spectral features 
are discussed in Sec. IIIIl An analytical expression for the 
optical conductivity <t{u) is derived in Sec. |TV] where we 
also investigate the different polaronic features and their 
dependence on the microscopic parameters. In Sec. IVlwe 
present a further simplified expression for (7{lu) valid in 
the lattice polaron regime. The main results are briefly 
summarized in Sec. IVII where also the consequences of 
including spin fluctuations (here neglected) is also dis- 
cussed. Finally, a detailed derivation of the analytical 
expression for the one-particle spectral function and the 
optical conductivity is reported in the Appendices. 



II. HOLSTEIN-i-,7 MODEL IN INFINITE 
DIMENSIONS 

In the following we consider the case of a single hole 
in an antiferromagnetic (AF) background interacting 
with local Holstein phonons. Using the linear spin- wave 



approximatiort 
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and neglecting terms that vanish 



write the Hamiltonian asi^ 

t /. +. . \ J(l - 22:) 



+ gY,hlh.,(b, + bl)+u;oY,blh 



(1) 



Here is the creation operator for boson spin defects, 
/i^ is the single spinless hole operator and x ~ (a^a) 
represents the density of spin defects which is finite at 
nonzero temperature. Note that in the thermodynamical 
limit the presence of a single hole does not affect the 
magnetic state, which can be thus evaluated (in the z ^ 
1 limit) in the absence of spin dynamics. The density 
of spin defects can be obtained from the magnetization 
m — 1 — 2x via the Curie- Weiss equation: 



m 



tanh 



(2) 



which defines a Neel temperature Tn = J/4. Concern- 
ing the electron-lattice interaction, we shall mainly fo- 
cus on the adiabatic regime loq ^ t which is relevant to 
the experimental systems of interest. In this regime, a 
dimensionless electron-phonon coupling constant can be 
defined as A = g"^ /uqI, the polaron energy in units of the 
hopping integral. 

An exact solution for the thermodynamical and the 
one-particle spectral properties of Eq. ([1]) at T = was 
obtained in Ref. [s^ in terms of a continued fraction. In 
order to derive the optical conductivity, the one-particle 
Green's function must be generalized to finite tempera- 
ture, which involves the following steps: i) one has to 
allow for thermally excited phonons; ii) the presence of 
thermally excited spin defects requires the introduction 
of a "spin-resolved" Green's function, to distinguish hole 
excitations created on sites with/without spin defects; 
iii) finally, the reduced magnetization introduces an ef- 
fective exchange coupling J = Jm < J. The details of a 
formal derivation of the one-particle propagator at finite 
temperature are reported in Appendix El we summarize 
here the main results. 

Following Ref. [13, we define Gifi{uj) — Gi{uj) as the 
Green's function for one hole created on a site in the 
absence of spin defects. A careful analysis (see Appendix 
[A]) shows that the Green's function Gi^i{uj) for a hole 
created on a site with a spin defect is simply Gi^i{tu) = 
Gi{u! + J). In addition, at finite temperature the Green's 
function Gi {uj) itself is defined as a thermal average over 
the phonons: 



Giu) 



It' ;V, 



-(3nujQ r^nn / 



- nwo). 



(3) 



in the limit of large coordination number z 3> 1, we can 



where (^""(tj) represents the propagation of one hole 
created on a site with n excited phonons and 2'ph is 
the single-site phonon partition function = 1/(1 — 
Q-P'^oy Following Ref. |38j, we can derive a self-consistent 
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expression for G"" in terms of a continued fraction. We where 
can write: 

I 



KM 



g-\LO-{n+l)uo) 



g-i(cj-(n + 2Vo)- 
"5^ 



(n + 3)5" 



(n-l)g^ 



(5) 



(6) 



{n - 2)g' 



Carrier motion and exchange interactions are taken into 
account by the bath propagator 

g-\Lo)^Lu~i:t{oj), (7) 

where the "hopping" self-energy is defined as 

Y^tiuj) = — - x)Gj{lu - J /2) + xGj{uj + J/2) . 

(8) 

FinaUy, the spin-defect averaged Green's function, which 
is the physicaUy probed quantity in photoemission exper- 
iments, is obtained as: 



G{uj) = (1 - x)G{lu) + xG{lo + J/2). 



(9) 



Note that the factors (1 — x) and x in front of Gjiuj— J /2) 
and Gj(io -I- J/2) account for the probability of a site to 
be, respectively, free or populated by a spin defect. Note 
also that the Green's function G(lS) appearing in Eqs. 
(I7l8p is a phonon averaged quantity, so that the solution 
of Eqs. involves the simultaneous self-consistency 

of all of the G"", which is much more computationally 
expensive than the solution at zero temperature. 

It is easy to check that in the absence of electron- 
phonon interaction Eq. ^ recovers the thermal Green's 
function for the pure t-J model as defined in Eq. (2.12) 
by Stumpf and Logani^ On the other hand, it should be 
stressed that the present solution, although it formally 
recovers the results for the pure Holstein model in the 
limit J ^ (paramagnetic case) [Eqs. (40)-(42) of Ref. 
[38| . is still described by a purely local Green's function 
Gi_j[ijj) = bi_jG(uj) due to the assumption of a classical 
(although disordered) spin background. The main draw- 
back of this assumption is thus that no coherent disper- 
sive peak is obtained in this framework even in the J — *■ 
limit of the t-J Holstein model, and, consequently, no 
Drude peak appears in the optical conductivity. Notwith- 
standing this limitation the present approach, as we will 



show below, is still quite able to reproduce in more than 
qualitative agreement the incoherent features of both the 
one-particle Green's function and of the optical conduc- 
tivity. 



III. ONE-HOLE SPECTRAL PROPERTIES 

Before discussing the optical conductivity, let us briefly 
present our results for the one-particle spectral function 
in the Holstein-i-J model at finite temperature. Pre- 
vious approaches have focused separately on thermal ef- 
fects either in the pure Holstein or in the pure t-J model. 
The temperature evolution of the hole spectral function 
piuS) = — (l/7r)ImG(a;) for the t-J model in the infinite 
dimensional limit has been analyzed in Ref. i27i . At 
T = 0, it consists of a series of (5-function magnon peaks 
whose distribution refiects the strength of the magnetic 
polaron: for sufficiently large J/t the profile is rapidly 
decaying with energy (reminiscent of the string picture 
of small magnetic polarons) whereas for small J/t it 
acquires a more symmetric shape, reducing to a semi- 
circular function in the limit J/t 0. The effect of a 
non-zero temperature within this context is to broaden 
each ^-function with a bandwidth Wx which is ruled by 
the thermal spin defect probability x. We shall term this 
effect the intrinsic magnetic broadening. Such intrinsic 
magnetic broadening Wx is however exponentially small 
for T/Tf^ ^ 1, and it becomes significant only close to 
Tn. a semi-circular shape is recovered also in the para- 
magnetic phase T > Tn, where J/t — 0. 

As the electron-lattice coupling is turned on, each 
magnon peak splits into several sub-peaks spaced by wq, 
reflecting the dressing of the hole by phononic excita- 
tions. In this context the strength of the electron-phonon 
coupling rules not only the number of phonon satellites 
but also its spectral weight profile. Just as in the pure 
Holstein model, while the number of phonon peaks is 
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FIG. 1: (color online) Temperature evolution of the total spec- 
tral function p{w) = — (l/7r)ImG(w) in the polaronic regime. 
Microscopical parameters: A = 0.7, uJo/t = 0.05, J/t = 0.4 
and T/Tn = 0.1, 0.5, 0.9 (corresponding to T/loq = 0.2, 1, 1.8). 
Numerical calculations have been done with a small imaginary 
frequency part rj = 0.007. 



At this point, it is useful to comment about the ef- 
fects of coherent hole propagation, that are implicitly 
neglected in our approach. These would induce a finite 
dispersion of order ~ J to the lowest energy magnon 
peakiii2£i2ii2^!^d£ It is clear that such dispersion would 
be visible only at sufficiently low temperatures and at 
moderate electron-phonon coupling strengths, when the 
energy scale J is smaller than both the intrinsic magnetic 
broadening and the Gaussian phonon spread, whereas in 
the opposite case it will presumably be hidden below a 
featureless background. These considerations give fur- 
ther support to the present DMFT approach in the pola- 
ronic and/or high temperature regime, where neglecting 
the coherent hole propagation would not affect signifi- 
cantly the spectral properties. As we shall see below, this 
is even more true for what regards the finite-frequency 
optical conductivity, where any dispersive peak would be 
convolved in any case with high-energy featureless struc- 
tures. 



IV. OPTICAL CONDUCTIVITY 



quite small in the weak coupling regime, in the polaronic 
state a large number of phonon satellites appear with 
a characteristic Gaussian profile. The envelope of the 
phononic fine structure has a spread which is governed 
by the energy associated with the lattice fluctuations: it 
is given by ^/Xujoi = g in the quantum limit, and in- 
creases VtxTt as T > wo/2.39 

In Fig. [T] we show the temperature evolution of a typ- 
ical polaronic spectral function, at A = 0.7, u!o/t — 0.05 
and J/t = 0.4. Throughout the paper, when not spec- 
ified, we shall take the hopping matrix element t as 
the energy unit. The temperatures considered here are 
T/Tn = 0.1, 0.5, 0.9, corresponding to T/lvq = 0.2, 1, 1.8. 
At the lowest temperature (T/Tn = 0.1) a phononic 
fine structure can be clearly seen, superimposed on the 
magnon peaks. The width of each phononic peak is due 
to the intrinsic magnetic broadening Wx described in Ref. 
[27I . It is exponentially small at this temperature, so that 
a small Lorentzian broadening rj = 0.007 has been intro- 
duced for clarity. On the other hand, the spread of the 
multiphonon structure gives rise to an overall width to 
the magnon peaks that is in good agreement with the 
expected value g — 0.19. 

Upon increasing the temperature to T /T^ — 0.5, two 
different effects are visible. First, the width Wx of each 
of the fine peaks increases due to the intrinsic magnetic 
broadening, leading to a much smoother curve (this ef- 
fect actually overcomes the small Lorentzian broaden- 
ing 77 introduced previously). Also, the overall spread 
of the multiphonon profiles increases due to the thermal 
phonon fluctuations, as expected for T ^ ujq/2. Finally, 
at T/Tn = 0.9, the system is so close to the paramagnetic 
phase that neither the phonon peaks nor the magnon 
structure can be resolved. 



The evaluation of the optical conductivity for a sin- 
gle hole is a delicate matter, which is only partly sim- 
plifled in the context of DMFT due to the absence of 
vertex corrections^^ A controlled procedure is derived 
in terms of an expansion of the inverse fugacity at fl- 
nite temperature, performing the limit /i — > —00 to en- 
force the thermodynamically vanishing particle density. 
We can thus deflne the optical conductivity per hole, 
cr(w) = lim„,^_,o ""(w; n/i)/n/i, which is a flnite quantity 
and which presents the same features as in the dilute (but 
finite) hole density limit. Applying this formalism one 
derives a similar expression as obtained in Refs. [T9ll27l . 
here adapted to take into account the electron-phonon 
interaction. Leaving once more the technical details in 
Appendix|Bl we report here the main results. The optical 
conductivity per hole is expressed as 



a{uj) 



4uj 



X xj5{uj + n + J/2) + {l-x)p{Lu + n-J/2) ,(10) 



where 



and 



--lmG(w), 

TT 



p{u)) 



J dcjc'^'^piuj)' 



(11) 



(12) 



The last line defines a "weighted spectral function", 
P^{lu), which represents thermally excited states and 
plays an important role in the determination of the opti- 
cal conductivity. 

Let us remark that, although Eqs. pT7)) - (fT^ are for- 
mally analogous to those for the pure t-J model, the 
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FIG. 2: (color online) Comparison between the optical con- 
ductivity obtained by our DMFT solution and Lanczos 
diagonalization in two dimensions on a finite cluster (Ref. Isol . 
arbitrarily scaled). A Gaussian broadening A — 3aJo/5 for 
ct((jj) has been employed in our DMFT analysis (see text for 
details). 
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electron-phonon interaction appears implicitly in them in 
the evaluation of the local Green's function G{oj). Note 
also that, in the paramagnetic hmit J ^ 0, Eqs. (fTU]) - 
(fT2| do not recover the results of the Holstein model, '^'^ 
because Eq. (|10p involves the convolution of two lo- 
cal rather than li.- dependent propagators. As discussed 
above, this is due to the classical treatment of the spin de- 
grees of freedom which does not allow for coherent trans- 
port, so that no Drude peak is recovered in the present 
analysis. This, however, has only a minor influence on 
the finite-frequency optical conductivity, which is domi- 
nated by local incoherent excitations. 

In order to assess the validity of the present treatment, 
we compare in Fig. the optical conductivity of the 
Holstein-t- J model in infinite dimensions, as described by 
Eqs. ifTUIl - P^ . with numerical calculations using Lanc- 
zos diagonalization for a single hole in the 2D Holstein-i- 
J model on a VTO x -\/T0 cluster.'^" For technical reasons 
(see discussion below), DMFT data are averaged with a 
Gaussian filter of amplitude A = 3a;o/5 such that phonon 
resonances are still well separated. Microscopic values are 
A = 1, J/t = 0.4, wo/t = 0.2, and T = (for Ref. 's^ and 
T — O.Oli ~ O.ITn for the present results. These values 
correspond to a case where the lattice/magnetic polaron 
is formed and incoherent contributions to the optical con- 
ductivity are indeed dominant. The good agreement of 
the overall shape confirms the feasibility of our analysis 
to investigate the finite frequency optical conductivity. 



A. Technical details 

As discussed in the previous Section, the one-hole spec- 
tral function consists of a set of narrow bands whose 
width Wx is controlled by the intrinsic magnetic broad- 
ening. This intrinsic width vanishes at low temperature 
due to the absence of spin-wave dispersion, making the 



nLiiJ I I I 1^ ' I I I ■ ^ 

1 2 3 1 2 3 

co/t co/t 

FIG. 3: (color online) Left panels: Gaussian averaged op- 
tical conductivity as a function of the sampling mesh Auj 
for different temperatures. Right panel: the corresponding 
Gaussian averaged optical conductivity in the presence of a 
disorder-induced bandwidth Wd / Alj = 40 as described by our 
scaling procedure. The microscopic parameters are A = 0.7, 
uo/t = 0.05 and J/t = 0.4. 



direct evaluation of Eqs. pUjl - P^ quite hard to per- 
form. Special care is thus needed in order to calculate 
numerically the optical conductivity. 

We shall be mainly interested in the features of the op- 
tical conductivity that are related to the electron-lattice 
coupling. We shall therefore retain the details of the spec- 
tra on the scale of the phonon frequency loq. From the 
practical point of view, the primary object of our analy- 
sis will be a Gaussian average of the optical conductivity 
with a Gaussian filter of amplitude A = 3a;o/54^ 
Note that, even though the Gaussian average preserves 
the total spectral weight of the original set of data, the 
accuracy of the final result will be limited by the finite 
frequency sampling Aw. In explicit terms, if the spac- 
ing Aw is larger than the intrinsic peak- width Wx^ the 
data sampling will probe the spectral features in a ran- 
dom way, yielding a highly inaccurate result for both 
the shape and spectral weight of the optical conductiv- 
ity. This is shown in the left panels of Fig. [3] where we 
plot the dependence of the (Gaussian averaged) optical 
conductivity on the sampling spacing Aw for different 
temperatures. At high temperature T/Tn = 0.7 the in- 
trinsic magnetic broadening Wx is large enough so that 
both the shape and the total spectral weight of the opti- 
cal conductivity are well captured even with a relatively 
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large mesh {Aui/t = 7.8 x 10~^). At lower temperature 
T/T^ — 0.5 however becomes so small that a much 
finer mesh (Aw/t = 1.9 x 10^'^) is needed in order to get 
accurate results. At T/Tn — 0.3, finally, no convergence 
is achieved even for the finest sampling mesh considered 
in this paper, Acj/t = 1.1 x 10^* (in Fig. [3l for graphical 
reasons, we plot curves only up to Aw/t = 1.9 x 10""^). 
Clearly, since the intrinsic peak-width Wx vanishes ex- 
ponentially at low T, the problem evidenced here cannot 
be solved by merely reducing the sampling mesh. 

To overcome this difficulty we add to the system a 
small uncorrelated disorder with semielliptic distribution 
of amplitude Wd,— which is able to yield a finite peak- 
width even in the zero temperature limit. We choose 
Wd — 40Aa; to assure a sufficiently dense mesh for an ac- 
curate sampling. We then scale Wd keeping fixed the 
ratio Wd/Aw = 40 in order to approach the correct phys- 
ical limit in the absence of disorder. Results obtained 
with a finite Wd are shown in the right panels of Fig. 
[31 No appreciable difference is visible at high tempera- 
ture T/Tn = 0.7 where the thermally driven magnetic 
broadening Wx is large enough to guarantee the conver- 
gence even in the absence of disorder. On the other hand 
our procedure provides a clear improvement already at 
T/Tn = 0.5 where the convergence as function of the 
sampling spacing is more easily achieved in the presence 
of disorder (note that the converged results coincide with 
the results for the most dense mesh in the absence of dis- 
order, showing that no spurious structures appear due 
to our scahng procedure). Finally, for T/Tn — 0.3 the 
disorder scaling procedure is the only way to guaran- 
tee convergence of both the shape and the total spectral 
weight of the optical conductivity. 



B. Results 

Using the above described procedure, we shall now 
concentrate on the evolution of the optical conductiv- 
ity in the adiabatic regime, fixing the ratio uJo/t — 0.05. 
This value is qualitatively representative of the cuprates, 
where the half-bandwidth i « 1.2 eV and typical optical 
phonon frequencies cjph ~ 60 meV. This regime is also 
the most interesting one from the theoretical point of 
view, since in this case the interplay between lattice and 
spin degrees of freedom has its most dramatic effects. 

Fig. [4] shows the evolution of the optical absorption 
at fixed J/t — 0.4 and low temperature T/Tm — 0.3, 
upon varying the electron-lattice coupling. At A = 0.2, 
the result is very reminiscent of the spectra calculated 
for the pure t-J model in Ref. 27. It consists of a se- 
ries of magnetic peaks, dominated by the sharp single- 
magnon peak located at w ~ J/2 and rapidly decay- 
ing at higher frequency. In this weak-coupling regime 
A — 0.2, the electron-phonon interaction simply gives rise 
to a multi-phonon fine-structure with a Gaussian pro- 
file. Each magnon peak acquires thus a phonon-driven 
width without modifying however the overall distribu- 
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FIG. 4: (color online) Optical conductivity versus A across 
the polaron crossover, at fixed J — 0.4, uJo/t = 0.05 and 
T/loo = 0.6 {T/Tn = 0.3). The two arrows at low and high 
energy mark respectively the first magnon peak at lj ~ J/2 
in the weak electron-phonon coupling limit A <C 1, and the 
position u) ~ 2Xt of the broad polaronic absorption expected 
at A 3> 1 (shown here for A = 1.5). 



tion of spectral weight. 

The main effect of increasing the electron-lattice cou- 
pling is a progressive shift of the spectral weight towards 
higher frequencies. This is an evidence of the forma- 
tion of the lattice polaron, which occurs through a grad- 
ual crossover in the presence of a finite loq. Increasing 
A also modifies the shape of the low-energy absorption 
edge, converting the sharp magnon-peak at A ^ into 
a smoother Gaussian lineshape, typical of polaronic ab- 
sorption. In the strong coupling regime (A ^ 1), char- 
acteristic of a small lattice polaron, the position of the 
maximum in the optical conductivity is expected to scale 
linearly as lu — 2Xt. This is in good qualitative agree- 
ment with our data reported in Fig. [3] where however 
the effects of a finite hopping integral t (Ref. 1341) and 
of the 1 /uj behavior at high frequency (see Ref. UM ) re- 
sult in a slight redshift of the maximum of the polaronic 
structure. 

The evolution of the optical conductivity with A can 
be understood by analyzing the building blocks of Eqs. 
(|llll2p . In Fig. EJa), we report both the spectral func- 
tion p and the weighted spectral function p" for two typi- 
cal values A = 0.2 and A — 0.7. In the first case (dark red 
curve), there is no energy separation between the single- 
hole excitations in p and the thermally excited states in 
p" . The low- frequency gap in the optical absorption seen 
in Fig. m arises due to the explicit shift of the spectral 
function by the quantity J/2 in Eq. (fTO)) . representing 
the energy cost to create one spin-defect as the hole hops 
in the AF background. It is now interesting to compare 
these features with the results for A — 0.7 [light green 
curve in Fig. [Slja)]. As we can see, the spectral func- 
tions p(w) for A = 0.7 and A — 0.2 are qualitatively simi- 
lar, the only major difference being the increased number 
of phonon satellites involved in each magnetic peak, re- 
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FIG. 5: (color online) Spectral function p{uj) and weighted 
spectral function p"{i^) for the different cases reported in the 
legend. For better readability, the weighted spectral function 
p" (uo) is reported on the negative axis. The dashed black line 
in panel (a) refers to the strong coupling (SC) approximate 
formula as discussed in Sec. [V] 



fleeting the increased number of phonons in the polaron 
cloud. In order to understand the modification of the op- 
tical conductivity we thus focus on the weighted spectral 
function . The latter undergoes a much more drastic 
change reflecting in an explicit way the signature of po- 
laron formation. Of particular relevance is the shift of 
to much more negative energies which characterizes the 
lattice trapping. Note also that the magnetic peaks that 
are clearly visible at A = 0.2 are completely washed out 
at A = 0.7, merging into a broad polaronic-like spectrum 
centered at higher "binding" energies. This change in 
the nature of the thermally excited states is at the origin 
of both the opening of a polaronic pseudogap and of the 
smoothening of the features observed in cr(a;). 

A qualitatively similar evolution is observed upon vary- 
ing J, which gives a clear illustration of the positive inter- 
play between magnetic and lattice polaron effects. This 
is shown in Fig. [S] where we report the optical conductiv- 
ity for different values of J at constant A = 0.7 and at the 
same temperature T /uoq — 0.6 as in Fig. 2] (note that this 
corresponds to different T/T-^ as Tn scales with J). Re- 
markably, even though A is kept constant, reducing the 
magnetic exchange from the initial value J = 0.4 leads 
to a gradual loss of the lattice polaronic features and to 
the undressing of the hole from multi-phonon excitations. 
This results in a shift of spectral weight from the pola- 
ronic peak to a magnon peak at lower frequency. The 
magnon peak is already visible as a shoulder in the data 




FIG. 6: (color online) Optical conductivity across the 
polaron crossover, as varying J for fixed A = 0.7, 
uJo/t = 0.05 and fixed temperature T/wo = 0.6. Note 
that this fixed temperature corresponds to different T/Tn 
for J/t = 0.4,0.3,0.2,0.17,0.15, respectively T/Tn = 
0.3,0.4,0.6,0.7,0.8. 



of Fig. [Hlat J/t = 0.3 and J/t = 0.2, and it emerges more 
clearly at lower values of J. At J — 0.17, in particular, 
both the magnetic peak and the broad polaronic band 
are visible in the absorption spectrum. The featureless 
nature of the absorption curves at low J/t can be as- 
cribed to the increasing disorder of the AF environment 
as J diminishes and T/Tn increases 

Once again, the comparison of p" and p for J/t = 
0.4 [light green curve in (S^a)] and for J/t = 0.15 [dark 
blue curve in [SJb)] allows to visualize in a simple way 
the loss of the polaronic features. Even in this case the 
most relevant quantity is the weighted spectral function 
whose main excitations, for J = 0.15, are shifted to 
much less negative energies than for J/t — 0.4, closing the 
gap between and p. Note in addition that, although 
less evident than at J/t — 0.4, a magnetic peak in the 
spectral function p is still visible even for J/t = 0.15. The 
convolution of p with p"^ gives rise to the small magnetic 
peak at w ~ J/2 observed in the optical conductivity. 

In Fig. [7] we report the evolution of the optical spectra 
versus temperature in the polaronic regime at A = 0.7, 
LUo/t = 0.05 and J/4 = 0.4. Although Fig. [7] can look 
at a flrst glance quite similar to Fig. 21 no shift of the 
peak of (j{u!) is observed here. Instead, there is a pro- 
gressive filling of the low frequency gap with temper- 
ature that is quite similar to what is observed in the 
pure Holstein model.^" This can be pointed out again 
by the comparative study of p with p^. For T/T-^ — 0.9 
[light brown curve in[5l^b)] the weighted spectral function 
p" [ijj) still presents an extended peak at negative energy 
Lo/t w —1.3, signalizing that the lattice polaron is not 
completely destroyed. However, the broadening of p"{u}) 
is now significantly enhanced as T ^ wq- This feature, 
along with a similar broadening of the spectral function 
p{ijj) due to the approaching of the paramagnetic limit 
T/T^ — > 1 (see Sec. IIIip . leads to a significant overlap of 
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FIG. 7: (color online) Temperature evolution of the optical 
conductivity for A = 0.7, J/t — 0.4 and uo/t = 0.05. 



the two spectral functions and to a continuous filling of 
the gap in the optical conductivity as the temperature is 
increased. 



V. STRONG COUPLING FORMULA 

In the previous section we have discussed the techni- 
cal difficulties involved in the calculation of the optical 
conductivity. Even adopting the proposed scaling proce- 
dure (cf. Sec. IIV A[) . the computational cost can still be 
quite demanding, especially in the strong coupling lat- 
tice polaron regime where the number of phonons to be 
taken into account scales as = Xt/uiQ. In this section 
we present a simple analytical formulation of the opti- 
cal conductivity which is valid precisely in the adiabatic 
(luq 0) small polaron regime and which involves almost 
no numerical effort. 

The basis of this simplified formulation is the polaronic 
nature of the weighted spectral function p" pointed out 
in Fig. [5] We have seen indeed that at moderate values 
of the electron-lattice coupling, just above the polaron 
crossover, the function p" becomes essentially indepen- 
dent of the exchange term J. A closer look shows that in 
the lattice polaron regime the spectral function can 
be well approximated with the result of the atomic limit 



1 



V2t 



: exp 



(uj + 2\tf 



2s2 



(13) 



where s = g/ ■\/tanh(cjo/2r). This function for the case 
A 0.7, J/t 0.4, uja/t = 0.05, T/Tn = 0.3 is shown in 
Fig. [5^ compared with the full numerical solution for the 
same parameters. The good agreement can be ascribed to 
the atomic nature of the thermally induced excitations, 
which do not hybridize with the hopping continuum de- 
scribed by the self-energy Eq. ^ (an equivalent result 
was already pointed out in Ref. |34| ). Eq. (|13p provides 
thus a simple analytical expression for p"{lo) which can 



be employed in Eq. (|TC)|) for the evaluation of the optical 
conductivity. 

A simple analytical approximation can also be derived 
for the one-hole spectral function p{lo) involved in Eq. 
([To)) . At this level, we are mainly interested in the overall 
shape of the optical conductivity cr(tj), disregarding the 
multi-phonon fine structure that will be anyway washed 
out at large A as shown in Fig. [J] In this perspective 
we can consider formally the limit ^ and use the 
approach devised in Refs. [38ll39l as well as in Ref. [2 
to treat the adiabatic limit. In particular, the sum over 
the discrete energy levels in Eq. ([3]) can be replaced by 
an integration over Gaussianly distributed local random 
energies v which account for the thermal/quantum fluc- 
tuations of the phonon field. The local propagator thus 
becomes^i^^'^ 



G{lo) 



dv 



(14) 



where P{iy) = (1/V2^)e"'' Z'^" is a Gaussian function 
with the same variance as in Eq. (|13l ) . Self-consistency 
is achieved by employing Eqs. ([7])- ([8]). The comparison 
between the approximate formula for p{lo) valid in the 
polaronic regime and the full numerical solution is also 
shown in Fig. [li for A = 0.7, J/t = 0.4, uJo/t = 0.05, 
T/T^ = 0.3. Note that, although the bath propaga- 
tor G{uj) does not depend explicitly on the local random 
energies u, it still depends parametrically on A and J 
through the self-consistency conditions Eqs. ©-(El), so 
that it still retains the relevant lattice and magnetic spec- 
tral structures. 

Once the approximate analytical expressions for the 
spectral functions p and p" are provided, we can now 
simply evaluate the optical conductivity by using Eq. 
(fTO|l . Let us stress that the numerical solution of the 
approximate expressions for p and p", defined in Eqs. 
(fT3)l - (fT4)) requires a remarkably smaller computational 
cost than the full numerical solution, especially in the 
lattice polaronic regime where Eqs. P^ - (|14p are valid 
and where the computational cost of the full numerical 
solution is highest. The comparison of the approximate 
formula with the full numerical solution as function of 
the electron-phonon coupling constant A is illustrated in 
Fig. [H The agreement is excellent in the small polaron 
regime A > 0.7 and it is still quite satisfactory even for 
moderate electron-phonon coupling A = 0.5. 



VI. SUMMARY AND CONCLUSIONS 

In this work, we have provided an analytical treat- 
ment for the optical conductivity <7{lu) of one hole in 
the Holstein-t-J model at finite temperature, in the 
limit of infinite dimensions. In this context a dynamical 
mean-field solution can be derived where the local self- 
consistent problem is solved exactly. Due to intrinsic lim- 
itations enforced in infinite dimensions, we are not able 
to account for the coherent propagation of a single hole in 
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FIG. 8: (color online) Comparison of the optical conductivity 
evaluated from the full numerical solution (light green lines) 
and by means of the strong coupling formula (dashed black 
lines). The agrement is accurate already at moderate values 
of A. 

the antiferromagnetic background. Nevertheless, we have 
shown how the incoherent optical processes, related to 
local excitations (Einstein phonons and local spin-flips) 
are well described in our approach, as confirmed by the 
comparison with numerical results obtained by Lanczos 
diagonalization. 

Our main aim has been to investigate the incoherent 
features of cr(aj) in an intermediate coupling region where 
the positive interplay between the magnetic and lattice 
degrees of freedom is more relevant, sustaining the for- 
mation of a spin-lattice polaron. In this context we have 
studied the evolution of polaronic features in the optical 
conductivity as function of the different microscopic pa- 
rameters, as the electron-phonon coupling A, the temper- 
ature T and the effective exchange energy J. We remind 
that J does not represent the bare exchange energy but 
rather a mean-field-like Weiss exchange coupling which 
depends on the local magnetization m. In the cuprates, 
for example, this parameter can be tuned by varying the 
hole doping starting from the parent AF phase. We have 
shown that the role of the electron-phonon coupling is 
twofold. On one hand it rules the formation of the lat- 
tice polaron, changing the incoherent part of cr(w) from 
a typical spin-polaron spectrum, characterized by mag- 
netic peaks and by an optical gap at J/2, to a broad 
lattice-polaron-like shape located at higher frequency, 
from which the magnetic peaks are essentially washed 
out. On the other hand it also tunes the amount of 
quantum lattice fluctuations, reflected in the emergence 
of multi-phonon satellite peaks. This flne structure sur- 
vives also when the lattice polaron is destroyed at small 
A and it gives rise to an effective broadening of the mag- 
netic peaks which can be much larger than the intrinsic 



broadening driven by the thermal magnetic fluctuations. 

The present approach allows us to distinguish between 
two different mechanisms leading to a suppression of the 
polaronic pseudogap: i) at intermediate values of the 
electron-phonon coupling, the reduction of the effective 
exchange energy J leads to a shift of spectral weight from 
high energy lattice polaronic features to low energy mag- 
netic excitations. This results in a closing of the pseudo- 
gap as the hole undresses from its lattice polaron cloud; 
a) Conversely, increasing the temperature within the po- 
laronic regime, gives rise to a filling of the pseudogap. 
Both mechanisms are actually observed in the optical 
spectra of the underdoped cuprates.^'* 

As a final remark, we briefly comment on the ro- 
bustness of our results in physical systems where quan- 
tum spin fluctuations are present, allowing for coherent 
motion of the holes. As discussed above, one of the 
main effects is the emergence of a dispersive pole with 
small spectral weight in the one-particle spectral func- 
tion. In the optical spectra, this gives rise to a Drude-like 
low-frequency response, but it does not affect the high- 
frequency incoherent part (moreover the coherent spec- 
tral weight is strongly reduced when the electron-lattice 
coupling is turned on). Regarding the high-frequency 
part, the intrinsic dispersion of the spin fluctuations is 
itself expected to smear the magnetic peaks. This would 
affect our results only in the weak electron-phonon cou- 
pling regime, where no lattice polaron is formed. In 
this case an intrinsic broadening of the magnetic peaks 
in the optical spectra due to the spin-fluctuation dis- 
persion should be considered for a quantitative analy- 
sis. On the other hand, the effects of the dispersion of 
spin-fluctuations are expected to be barely visible in the 
lattice polaron regime, where the smearing due to the 
multi-phonon satellite structure around each magnetic 
peak dominates. 
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APPENDIX A: FINITE TEMPERATURE 
GREEN'S FUNCTION OF THE HOLSTEIN-t-J 
MODEL IN INFINITE DIMENSIONS 

In this Appendix we provide a detailed derivation of 
the Green's function of one hole in the Holstein-t- J model 
at finite temperature in the infinite dimensional limit. A 
formal derivation for the pure t-J model was discussed 
in Ref. |23, while the derivation for the full Holstein-i- 
J model at T = was provided in Ref. [H. On this 
ground, here we limit ourselves to the derivation of the 
finite temperature self-consistent equations in terms of 
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continued fractions. ping term while Hi^ contains all the other, purely local, 

Let us start by writing the Hamiltonian in Eq. ^ as contributions. We also define the Green's function as: 
H = Ht + Hi^, where Ht represents the non-local hop- 

I 



it) = -*^(i)^T^ E Mj#.; s., n.| e-'^^e-^^'cle^^'c. |n„ = 0; {n},^., {s},^.) , (Al) 



{n}.,{4. 



r 



where \ni,Si\{n}j^i,{s}j^i) denotes the state with rii 
phonons and Si spin defects (si = 0: no spin defect, 
Sj = 1: spin defect) on the site i and with a generic set 
{nlj^i of phonons and {s\j^i of spin defects on all the 
other sites. With these notations Ci — h\ when si — 
and Cj — h\a.i when Si — 1. In addition, Z'^{N) is the 
total partition function of the system in the absence of 
holes, TV being the number of sites; in this case Ht does 
not contribute and Z^{N) factorizes into a phonon and 



J 



a spin part as Z'^{N) = Zp'^(7V)Z"p'"(7V). Each of them 
can be in addition factorized with respect to the site in- 
dex, e.g. Z^P'"(7V) = ^r"nj^,^7'". Similar consid- 
erations hold true for the exponential terms e~^^e~'^* 
which apply on the states with no holes on the right side 
of Eq. (jAip . Reminding Hij\ni,Si) — riiUJo + SiJ/2, we 
can write, after few straightforward steps, in the Fourier 
space: 



G.(^) = 



-f3[muo^s^J/2] 
^ph ^spiii 



E 

{n}j^i,{s}j^i 



-f^J2j^il^j'^0+SjJ/2] 



X {{s}j-ii,{n}j-ii;s„n^\cl ~— -c, \n„ s,; {n}j^i, {s}j^i) 

u! + riiUjQ + SiJ / Z — H 



This can be rewritten as 



(A2) 



(A3) 



where 



E 



-l3siJ/2 



^spm 



E 



a-/3Z;j^i["3'^0 + Sj.//2] 



yph r^spiii 



uj + Si J/2 — H 



Ci Si]{n}j^i, {s}j^i) . 



(A4) 



Similarly, the factor ps — e PsiJ /2 j ^spm ^^gj^gg j-j^g local population of spin defects which can be evaluated within 
the mean-field theory enforced by the infinite dimensional limit, namely pi — x, po = 1 — x. We have thus: 



crM = (1 - x)Gi:sH + ^G™i'H, 



with 



Gzs^) = E 



-PJLj^i [njU>o+Sj J /2] 



vph /T-spin 



{{s}j^i,{n}j^,]Si,n,\ 



1 



Finally, using the definition of q, one can see that G""(w) = G™{u) 4- J/2), and we can write 



(A5) 



ni, Si; {n}j^i, {s}j^i) (A6) 



or (to) = (1 - x)Gr(^) + oocrito + Jm, 



(A7) 
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where 

g-/3E3^J"j'^o+Sj J/2] 2 

= 22 ^ph ^spin 7T|n,;{n}j^„{s}j^,). (A8) 

f-ir-i Zj ■ I ■ Zj ■ I ■ ^ ^ 



Here {n}j^i, {s}j^i) denotes the state with rij 
phonons on site i, {s}j^i being the phonon/spin 

configurations on the sites j 7^ i with one /lo/e on the site 



We can write Eq. (|A8p using the short-hand notation: 

1 , , 



Pi (.rii 



(A9) 



uo — H 

where \ni) = jn^; {n\j^i, {s}j^i), and where the operator 



Pi (n. 



E 



(AlO) 



denotes the average over the spin and phonon configura- 
tions on aU the sites but i. The operator Pi represents 
a direct generahzation of the quantity P{s) introduced 
by Stumpf and Logan in Ref. [23 to include the phonon 
degrees of freedom. It is also convenient to note that (for 



Pi (n-i 



-/3[njCJo+Sj J/2] 



ph yspiii 



zrz. 



X (5 j , Tij , 71^ I ... 1 72^ , Tij 1 Sj) . 



(All) 



where Pij is defined in similar way as Pi as the average 
over the spin and phonon configurations on all the sites 
except i and j. 

Having introduced the necessary definitions, from now 
on we can follow the derivation in Ref. |32 properly 
adapted to the finite temperature case. In particular we 
can introduce the local Green's function g""(a;) defined 
as the atomic t — limit of Eq. (jA9p . Note that in the 
atomic limit A = 1 so that 



1 



LU — nujo 



(A12) 



We can also generalize Eq. (|A9[) for off-diagonal local 
phonon matrix elements: 



1 



\Pi 



(A13) 



whose analytical expression will be provided In Eq. 

ISM- 

We can now employ the standard relation 



1 



1 



1 



LJ- H 



1 



1 



Ht 



1 



1 



.{MA) 



which, on a classical spin background, gives rise to the 
retraceable path approximation. We have thus 



Y,-g7{u)tf^{u)Grmi^) 



where sj^^'(w) represents the dynamics of the hole af- 
ter hopping on the neighboring site j. Since, in the 
leading term of a 1/d expansion, the hopping process 
and the further dynamics of the hole do not involve 
the phonon degrees of freedom on site i, such evolu- 
tion occurs in the presence of p phonons on the site i, 
which are reflected in an shift of the frequency argument, 
S^-^'''(tj) = Sj(ijj — ploq). Note that Sj(ijj) still contains 
the full thermal average on the phonons at site j, so that 
it is related to the thermally averaged Green's function 
defined in Eq. (jA3p . In addition, T,j{ui) will depend also 
on the initial spin configuration Sj at site j: hopping to 
a site j free of spin defects will create a spin defect on 
the site i, while hopping to a site i with a spin defect will 
restore the initial magnetic background at site i destroy- 
ing a spin defect at i. Using these considerations we have 
thus: 
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^p.,G,(a> + (2s,-l)j72) 



= — ^{1-x)Gj{lu~J/2) + xGj {uj + J/2)J (A16) 

and we can write Eq. (jAlSp in the compact form: 

G""(c.) = r^C^) -Y.~97i^)M^-P^o)&"^{iJiM7) 
p 

where we have dropped the unnecessary site indices and 
we have added the index to denote a self-energy term 
arising from the hopping processes. 

We can indeed employ now the usual procedure to ob- 
tain the electron-phonon self-energy in terms of a contin- 
ued fraction. In particular, in addition to the diagonal 
elements of Eq. (|A12p , we specify also the non-diagonal 
ones which read: 



= [lU - nUJo] Sn,m + gX^ 



(A18) 



Using the standard derivation, we can thus write the 
Green's function G™{uj) as 



G""(tj) 



1 



g~^{oj~nu;o)~K\{oj)-^:^^Juj) 
where Q~^{ll!) = uj ~ St(w), and 



, (A19) 
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KM , , 2 > (A20) 



f;-i(^-(n + lH) ^ ^ / ^ow 

g-^(c.-(n + 2H)- ^"^^^^ 



and 

r,2 



K^.M = ■ (A21) 

g [u- nujo) 



r-it I 1^ N ("--2)5 



APPENDIX B: OPTICAL CONDUCTIVITY 



In the previous Appendix we have derived an exact 
expression for the Green's function of a single hole in 
the Holstein-i-J model in infinite dimensions. Here we 
investigate within the same framework the optical con- 
ductivity per hole, cr(a;), in the zero density limit. To this 
aim we provide an alternative derivation with respect to 
Ref. 11^. The main advantage of the present approach is 
to deal at the same level with both the spin and phonon 
degrees of freedom, allowing thus for an immediate gen- 
eralization of the t-J model to the Holstein-t- J model. 
As a result we obtain a final expression of the optical 
conductivity as a functional of the one- hole Green's func- 
tion which is formally similar to the one of Ref. [23 but 
where the local Green's function takes into account the 
electron-phonon interaction. 

The formal way to derive the optical conductivity (t{uj) 
for a single charge, is to consider the limit a{u!) = 
liTHn^^o a{Lo;nh)/nh, where a{uj;nh), Uh are quanti- 
ties defined in the grand-canonical ensemble (O) — 
Y^Nn^'^^^'^^^^O^Nn /^G.c, and where Nh it the total 
number of charges. The limit n/j — > is enforced by ex- 
panding to lowest order in terms of the fugacity z = e^^^ 
(/i — !■ —oo). In this limit only the subscctor N]-^ = 1 
survives both in (j{Lj]nh) and in rih- 

Let us first consider the hole density, which we can 
write as: 



nh 



h]h, 



a 



(Bl) 



where |a) is a complete set of eigenstates with eigenvalues 
Ea of one single hole (subspace Nh = 1). Here, since the 
state I a) must contain one hole at site i, the hole- number 
operator is simply defined as Nh.i = h\hi without any 
spin defect. Inserting now into Eq. (|Bip a complete set 
of eigenstates I7) in the subspace A'';, = (no hole) and 



a (5- function, we obtain 

Uh = e^'^ I duje-'^'^ 



X 



■i,a,7 



-PE-, 



B2) 



Note that only the I7) which do not contain any spin 
defect at site i would contribute to Eq. (IB2|1 . otherwise, 
after the hole creation, we would end up with a state 
with one hole and one spin defect present, which is for- 
bidden in the Hilbert space. We have thus: 



Uh = Ne^^'poj dwe"^'^ 

-HE^ 



S{uj-Ea+E^) 



where po is the statistical probability to have a site 
with no spin defect and where now only the I7) with 
no spin defect at the site i are selected. The square 
bracket in Eq. (|B3[) is just the local spectral function 
p{uj) — — (l/7r)ImG'o(i^), namely 



a, 7 



-I3E-, 



7 



d{Lu-Ec, + E-,), (B4) 



[remind that Go{uj) — G{uj)], so that we have 



(B5) 



Let us turn now to the optical conductivity or, 
more precisely, to the current-current response function 
Tl{i-u) which is related to cr(uj; Uh) through the relation 
a{uj;nh) = — ImH(w -I- iS)/uj. According to the previ- 
ous argumentation, in the limit Uh — > we can limit our 
analysis to the Nh = 1 subspace and write: 



(B6) 
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where r is the imaginary time in the Matsubara space 
and J is the current operator to be defined below. 

After usual manipulations we can write in the Fourier 
space: 



\ - e 



J-. 



H + E, 



-J 



(B7) 

where ujm = 27rmT are bosonic frequencies and where we 
remind that \a) are eigenstates in the = 1 subspace. 
The current density operator can be written as = 

{itjl-Jz) 'Ylii^i j) {^c\cj — c]ci^ , where we are summing ex- 
plicitly on all possible directions (this prescription com- 
pensates for the well-known vanishing of the current- 
current response in infinite dimensions'^). As discussed 
in Refs. [s^fssl and as appearing in Eq. ([I]), cjcj = hihjGj 
if there is a spin defect on the site j {sj = 1), while 



cjcj = alhih^j if there is no spin defect on the site j 
{sj = 0). Due to the classical nature of the magnetic 
background, it is easy to realize that a retraceable path 
approximation is enforced also in the current-current re- 
sponse function just as in the one-particle Green's func- 
tion. We obtain thus: 



n(icj,„) 



^4^^ 



- H 



As usual, we can now insert twice in Eq. (|B8[) the 
identity operator |7)(7|, where I7) are eigenstates in 
the subspace without any hole. We have thus: 



n(iw„) 



4z 



E ■ 

a, 7, 7' 



-/3-Eq 



P{H-E^) 








C," 


-'') 


(y 





(B9) 



Let us summarize the physical meaning of Eq. (jB9p . The eigenstate \a) contains one hole at the site j, and Eq. 
(|B9p describes the hopping of the hole to site i. If no spin defect is initially present at i {si = 0), this process involves 
the creation of a spin defect at j, and hence Ci = h^, = a^hj. On the other hand, in the alternative case where a 
spin defect is initially present at i {si = 1), the hopping process destroys the spin defect at the site i, so that Ci = aih\ 
and c] = hj. 

Let us consider for the moment the first case. We have thus: 



1X0(1 w„ 



4z 



a, 7, 7' 



E 



1) { 1 



-^I3(H~E^) 



-H + E„ 



alh 



(BIO) 



Let us now ask ourselves the following question: how 
much does the one-hole state \a) differ from the "free- 
like" state I7) in the absence of holes? It is clear that 
in the t — case only the phonon-spin configuration on 
the site j is affected by the presence of the hole, whereas 
all the other sites would be unaffected. In the presence 
of hole dynamics {t ^ 0), however, all the other sites 
are in principle affected. We remind also that, because 
of the classical magnetic background, the hole dynam- 
ics obeys a retraceable path approximation just as in a 
Bethe lattice, as depicted in Fig. [9l Let us consider now 
a given specific link of two nearest neighbor sites, 
In a true Bethe lattice, the whole system can be divided 
in two subspaces, (1) and (2), connected by the hopping 
term tij oc tj'Jz. This means that, to leading order, 
the presence of one hole on the site j would not affect 



the subspace (1). Similar considerations hold true for a 
generic lattice under the retraceable path conditions: ad- 
ditional links between the two subspaces (dashed lines in 
Fig. [9]), will contribute only to 0{l/y/z) and they can 
be neglected in the z — > 00 limit. For practical purposes 
we can thus split the total Hamiltonian as H — Hi + H2, 
where Hi accounts for the phonon-spin degrees of free- 
dom of the subspace (1) (not including j), while H2 con- 
tains the phonon-spin degrees of freedom of the subspace 
(2) (not including i). In a similar way the eigenstates 
\a), I7) can be written (in the leading order of a 1/z 
expansion) as \a) — \oli) ® \ol-2) and I7) = I71) ® I72). 
Employing these results, and noting that (Q\h},^aj\'-f) = 

(ai|/i]aj|7i)fo2,72, {i\a\hj\a) = (7;|a]/ij|ai)57^,a2, 
(Cf\hi . . . h\\i) = (j2\hi ■ ■ ■ /j1|72>'57i,7i: we have thus: 
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FIG. 9; Schematic picture of the hole dynamics in the retrace- 
able path approximation for z ^ oo. Given a hole at site 
j [subspace (2)], the probabihty to affect the subspace (1), 
through the Unk or other links (dashed connections), is 
0{^l \fz). In the leading order 1/z — » the subspace i/i, Hi 
are thus independent. 



4^ 



E 

0'>»,a;i,7i 



ZJ 



ai 



h]aj 



71 



^2 



a.2 



I _ Q-l3{E-n+H2-Eai-Ea2) ^ 



a2 



(Bll) 



where Z]', Zj are the partition functions of the corresponding subspaces and J2{j)i denotes a sum over the z nearest 
neighbors of the site i. Performing the analytical continuation ioJm u) + i5, and introducing once more appropriate 
^-functions, we end up thus with: 



CTo(w;n/j) 



Imno(a; + i6) 



tVf^ [1 - e-^'^] 



. 0'>i,ai>7i 



0L\ 



h]aj 



71 



E 



{a2\hi\X2)\^SiLo + i}-Ex,+E^,) 



(B12) 



where IA2) are cigcnstatcs in the Nh = 1 subspace, po = 
1 — X, pi = X are the statistical probabilities to have no 
spin defect and one spin defect, respectively, and where 



Pi 



M = E' 



a,7 



-f3E^ 



a 



h\ai 



7 



5{uj -Ea + EjIBU) 



The same derivation can be now employed for the case 



when a spin defect is present on the site i. After few 
straightforward calculations, we obtain: 



ai{uj;nh) = -PoPi 

X y dQe-^^po{fl)pi{iJ + fl). (B14) 
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Summing the two contributions (|B12p . (|B14p . and re- 
minding po(<^) = p{^), = p(w+J/2), after a change 
of variable we obtain: 

cr(w;n,i) = o-o(w;n,i) +ai{uj;nh) 
t^TTNe^f" [1 - e-P'^] po 



X pip{lj + n + ,1/2) + pap{uj + n- J 12) (pi5) 



where we made use also of the relation pie^'^/'^ = po. 
Finally, dividing Eq. (|BT5)) by dBSj we obtain the di- 
mensionless quantity Eq. (jlOp. Note that, dividing Eq. 
(|B15p by (|B5p . the common factors A^e''^ cancel out, so 
that the limit jj, —oo is well defined. 
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